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ABSTRACT 

The flow of energy through the solar atmosphere and the heating of the Sun's outer regions are 
still not understood. Here, we report the detection of oscillatory phenomena associated with a large 
bright-point group that is 430,000 square kilometers in area and located near the solar disk center. 
Wavelet analysis reveals full- width half-maximum oscillations with periodicities ranging from 126 to 
700 seconds originating above the bright point and significance levels exceeding 99%. These oscilla- 
tions, 2.6 kilometers per second in amplitude, are coupled with chromospheric line-of-sight Doppler 
velocities with an average blue shift of 23 kilometers per second. A lack of cospatial intensity oscil- 
lations and transversal displacements rules out the presence of magneto-acoustic wave modes. The 
oscillations are a signature of Alfven waves produced by a torsional twist of ±22 degrees. A phase shift 
of 180 degrees across the diameter of the bright point suggests that these torsional Alfven oscillations 
are induced globally throughout the entire brightening. The energy flux associated with this wave 
mode is sufficient to heat the solar corona. 
Subject headings: Published in Science 20-March-2009 



1. REPORT 

Solar observations from both ground-based and space- 
borne facilities show that a wide range of magneto- 
acoustic waves (1,2) propagate throughout the solar at- 
mosphere. However, the energy they carry to the outer 
solar atmosphere is not sufficient to heat it (3). Alfven 
waves (pure magnetic waves), which are incompressible 
and can penetrate through the stratified solar atmo- 
sphere without being reflected (4), are the most promis- 
ing wave mechanism to explain the heating of the Sun's 
outer regions. 

However, it has been suggested that their previous 
detection in the solar corona (5) and upper chromo- 
sphere (6) is inconsistent with magnetohydrodynamic 
(MHD) wave theory (7, 8). These observations are best 
interpreted as a guided-kink magneto-acoustic mode, 
whereby the observational signatures are usually sway- 
ing, transversal, periodic motions of the magnetic flux 
tubes (7, 9). Numerical simulations (10) show that 
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subsurface acoustic drivers and fast magneto-sonic kink 
waves (11, 12) can convert energy into upwardly prop- 
agating Alfven waves, which are emitted from the so- 
lar surface. These numerical simulations are also in 
agreement with current analytical studies. In particu- 
lar, it has been shown that footpoint motions in an axi- 
ally symmetric system can excite torsional Alfven waves 

(13) . Other Alfven wave modes may exist, although 
these are normally coupled to magneto-sonic MHD waves 

(14) . In the solar atmosphere, magnetic field lines clump 
into tight bundles, forming flux tubes. Alfven waves 
in flux tubes could manifest as torsional oscillations (7) 
that create simultaneous blue and red shifts, leading to 
the non- thermal broadening of any isolated line profile, 
and should thus be observed as full- width half-maximum 
(FWHM) oscillations (15). A promising location for the 
detection of Alfven waves is in the lower solar atmo- 
sphere, where they can be generated by the overshooting 
of convective motions in the photosphere (16). Here, we 
report the detection of substantially blue-shifted plasma 
and FWHM oscillations originating in a large conglom- 
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Fig. 1. — Simultaneous images in the (left) Ha continuum (photosphere) and (right) Ha core (chromosphere) obtained with the SST. 
The conglomeration of bright points within the region we investigated is denoted by a square of dashed lines. The scale is in heliocentric 
coordinates where 1 arc sec ~ 725 km. 



eration of magnetic bright points. 

We used the Swedish Solar Telescope (SST) to image 
a 68-by-68-arc sec region on the solar surface positioned 
near the disk center. Using the Solar Optical Universal 
Polarimeter (SOUP) (17) and high-order adaptive optics 
(18), we obtained narrowband images across the Ha ab- 
sorption profile centered at 6562.8 A. We observed with a 
cadence of 0.03 s to obtain 89 min of uninterrupted data. 
Because SOUP is tunable, we sampled the complete Ha 
line profile using seven discrete steps. The wavelength 
intervals we chose became increasingly narrow toward 
the line core in order to enable an accurate determina- 
tion of the line characteristics, such as Doppler velocities, 
FWHMs, and intensities. Our images have a sampling of 
0.068 arc sec per pixel, which corresponds to ^110-km 
resolution (two pixels) on the solar surface. 

By using the multi-object multi-frame blind deconvo- 
lution (MOMFBD) (19) image restoration technique to 
remove the small-scale atmospheric distortions present 
in the data, we achieved an effective cadence of 63 s for 
a full line profile. We acquired 85 complete scans across 
the Ha line profile in addition to 595 simultaneous con- 
tinuum images. For each of the 85 profile scans, every 
pixel of the 1024-by- 1024-pixel 2 charge-coupled device 
contains information acquired at a particular wavelength 
position. Therefore, we obtained a total of 8.9 x 10 7 indi- 
vidual Ha absorption profiles, covering the full 68-by-68- 
arc sec field of view, during the 89-min duration of the 
data set. We fit a Gaussian distribution to each of the 
observed Ha profiles to obtain values for the integrated 
intensity and FWHM. To determine the line-of-sight ve- 
locity, we compared each measured central wavelength 
position with the rest-frame Ha profile core at 6562.8 A. 
We created time series for intensity, line-of-sight velocity, 
FWHM, and wavelength-integrated data cubes and used 
fast Fourier transform and wavelet routines to analyze 
them. 

The SST field of view shows a range of features, includ- 
ing pores, exploding granules, and a multitude of bright 
points (Fig. 1), a large conglomeration of which is lo- 
cated at heliocentric coordinates (—10 arc sec, 10 arc sec) 



or N07E01 in the solar north-south-east- west coordinate 
system. We selected a 10-by-10-arc sec box surrounding 
the bright-point group (BPG), which occupies an area of 
430,000 km 2 , for further investigation. The line-of-sight 
Doppler velocities associated with this BPG show blue 
shifts with an average value of 23 km s _1 . There is no 
evidence of periodic trends in either intensity or line-of- 
sight velocity; the intensity of the BPG is constant, with 
minimal variation during its 53-min lifetime. 

Wavelet analysis shows that FWHM oscillations with 
significance levels exceeding 99% occur within the spa- 
tially averaged BPG (Fig. 2). We detected FWHM oscil- 
lations as low as the Nyquist period (126 s) throughout 
the duration of the data set, with the strongest detected 
power originating in the 400-to-500-s interval. These os- 
cillations are located directly above the large BPG, en- 
compassing a near-circular shape that is cospatial with 
the detected Doppler velocities, and are apparent in all 
FWHM time series. This shows that powerful coherent 
periodicities are present throughout the surface of the 
BPG. We detected oscillations until the BPG fragmented 
into a series of smaller bright points after 3150 s. 

Numerical simulations based on three-dimensional 
magnetoconvection show that the bright points that were 
observed in the wing of the Ha line profile correspond 
to magnetic field concentrations measured in kilogauss 
in the photosphere (20). The canopy structure seen in 
the Ha core images reveals a wealth of flux-tube struc- 
tures, with many securing anchor positions in the pho- 
tosphere directly above the BPG (Fig. 1). The coinci- 
dence of bright-point structures with high magnetic field 
concentrations implies that MHD waves are likely to be 
present (21). However, the chromospheric brightening is 
of much larger physical size than the underlying photo- 
spheric BPG. Because the observations were made very 
close to the center of the solar disk, an increase in phys- 
ical size between the photosphere and the chromosphere 
can be interpreted as an expansion of the photospheric 
flux-tube bundle as a function of atmospheric height (22). 
A comparison of the maximum diameter of the bright 
point at each height in the atmosphere suggests an ex- 
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Fig. 2. — A wavelength- versus-time plot of the Ha profile show- 
ing the variation of line width at FWHM as a function of time. 
The arrows indicate the positions of maximum amplitude of a 
420-s periodicity associated with the bright-point group located 
at ( — 10 arc sec, 10 arc sec) in Fig. 1. The torsional motion of the 
Alfvenic perturbations creates nonthermal broadening that is visi- 
ble in the Ha line profile. The peak-to-peak velocity is ~3.0 km s _1 
(?^65 mA). For an inclination angle of 35°, the absolute velocity 
amplitude is ~2.6 km s -1 . 

pansion of ^1300 km; a height separation of ^1000 km 
and a symmetric expansion around the bright-point cen- 
ter suggest a flux-tube expansion angle of ~33°. Addi- 
tionally, an offset of ~700 km between the centres of the 
BPG at photospheric and chromospheric heights suggests 
a magnetic flux-tube tilt angle of ~35° from the vertical. 

Alfvenic fluxes are predicted to be at their strongest 
in the regime of high magnetic field strength and mod- 
erately inclined waveguides (10). Because of their in- 
compressibility, they exhibit no periodic intensity per- 
turbations. Thus, the observational signature of a tor- 
sional Alfven wave propagating with a velocity compo- 
nent along the observer's line of sight will arise from its 
torsional velocities on small spatial scales (8). These tor- 
sional velocities are responsible for the FWHM oscilla- 
tions we observed (Fig. 3). The line-of-sight velocity am- 
plitude of wl.5 km s _1 and the inclination angle of ~35° 
indicate an absolute Alfvenic perturbation amplitude of 
^2.6 km s _1 . Because the circumference of the photo- 
spheric bright point [where torsional Alfven waves are 
believed to be generated (16)] is on the order of 2800 km 
(55 pixels), a torsional twist of ±22° is sufficient to gen- 
erate the observed wave motion. 

The moderate inclination angle of the flux tubes, cou- 
pled with the detection of substantially blue-shifted ma- 
terial and strong FWHM oscillations, is evidence of the 
presence of torsional Alfven waves. For a typical pho- 
tospheric internal waveguide electron density (23) of 
n\ ~ 10 16 cm -3 and a magnetic field strength (20) of 
1000 G, the Alfven speed within a cylindrical flux tube is 
estimated (14) to be ~22 km s _1 . This value is above the 
speed of sound in the upper photosphere/ lower chromo- 
sphere (24) and is consistent with the blue-shift velocity 
we determined. 

We took a slice through the center of the bright point 
and analyzed the stability of opposite edges of the BPG 
as a function of time. This was performed by examining 
any displacements of the BPG from its initial position 
at the start of the observing sequence (fig. SI). The 
bright point moves less than one pixel during the first 
3150 s. As the BPG begins to fragment after 3150 s, the 




Fig. 3. — Expanding magnetic flux tube sandwiched between 
photospheric and chromospheric intensity images obtained with the 
SST, undergoing a torsional Alfvenic perturbation and generating 
a wave that propagates longitudinally in the vertical direction. At 
a given position along the flux tube, the Alfvenic displacements are 
torsional oscillations that remain perpendicular to the direction of 
propagation and magnetic field outlining constant magnetic sur- 
faces. The largest FWHM will be produced when the torsional 
velocity is at its maximum (at zero displacement from the equilib- 
rium position). The figure is not to scale. 

motions of the bright-point edges increase substantially. 
However, we did not find periodic motions of the bright 
point, particularly during the initial 3150 s when the 
FWHM oscillations were detected. A magneto-acoustic 
wave mode would produce observable periodicities in in- 
tensity, similar to those caused by the periodic contrac- 
tions when viewed along the flux tube, that are asso- 
ciated with sausage- mode waves (25). A sausage- mode 
wave is caused by the axially symmetric expansion and 
contraction of magnetic flux tubes (14). Kink- mode os- 
cillations are generated through a bulk motion, whereby 
the whole flux tube is displaced from its original position 
in a periodic fashion. Therefore, magneto-acoustic waves 
cannot explain our observations. 

A torsional Alfvenic perturbation should produce a 
FWHM oscillation that is 180° out of phase at oppo- 
site boundaries of the waveguide (8). The relative time- 
averaged oscillatory phase as a function of distance across 
the ^2200-km diameter of the bright point shows that 
opposite sides of the bright point display oscillatory phe- 
nomena that are indeed 180° out of phase (Fig. 4). This 
is consistent with current torsional Alfvenic wave models 

(26) . 

We estimated the energy flux of the observed waves 
using E = pv 2 VA, where p is the mass density of the 
flux-tube, v is the observed velocity amplitude, and 
va is the Alfven speed (6). For a mass density of 
p w 1 x 10 ~ 6 kg m -3 , derived from a quiet- Sun chromo- 
spheric model (23), the energy flux in the chromosphere 
is E « 15000 W m -2 . At any one time, it is estimated 

(27) that at least 1.6% of the solar surface is covered by 
BPGs similar to that presented here. Thus, combining 
the energy carried by similar BPGs over the entire solar 
surface produces a global average of 240 W m -2 . Alfven 
waves with an energy flux of ^100 W m -2 are believed 
to be vigorous enough to heat the localized corona or to 
launch the solar wind when their energy is thermalized 
(6, 28). Therefore, a transmission coefficient of ~42% 
through the thin transition region will provide sufficient 
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Fig. 4. — Average phase difference of FWHM oscillations plotted 
as a function of distance across the diameter of the bright point. 
The triangles denote the locations where a measurement is made, 
and the error bars indicate the phase variance in the temporal do- 
main. The phase at km is used as a reference, with all other 
phases plotted relative to this value. Spatial coherence of FWHM 
oscillations across the BPG ranges between 90 and 100%, suggest- 
ing that the BPG is acting as a coherent waveguide. The phase 
difference increases around the midpoint of the BPG with opposite 
sides of the waveguide, indicating out-of-phase oscillatory phenom- 
ena, which is as predicted for a torsional Alfvenic perturbation. 

energy to heat the entire corona. Regions containing 
highly magnetic structures, such as bright points, should 
possess even higher mass densities (29). In this regime, 
the energy flux available to heat the corona will be sub- 
stantially higher than the minimum value required to 
sustain localized heating. 
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